arXiv: 1503.07267v3 [cond-matsupr-con] 25 Dec 2015 


Fractional Flux Plateau in Magnetization Curve of Multicomponent Superconductor Loop 


Zhao Huang and Xiao Hu^’^ 

^ International Center for Materials Nanoarchitectonics (WPI-MANA), 

National Institute for Materials Science, Tsukuba 305-0044, Japan 
^Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571, Japan 

(Dated: December 29, 2015) 

Time-reversal symmetry (TRS) may be broken in superconductors with three or more condensates interacting 
repulsively, yielding two degenerate states specified by chirality of gap functions. We consider a loop of such 
superconductor with two halves occupied by the two states with opposite chiralities. Fractional flux plateaus 
are found in magnetization curve associated with free-energy minima, where the two domain walls between 
the two halves of loop accommodate different inter-component phase kinks leading to finite winding numbers 
around the loop only in a part of all condensates. Fractional flux plateaus form pairs related by the flux quantum 
Oo = hcjle, although they individually take arbitrary values depending on material parameters and temperature. 
This phenomenon is a clear evidence of TRS broken superconductivity, and in a general point of view it pro¬ 
vides a novel chance to explore relative phase difference, phase kink and soliton in ubiquitous multi-component 
superconductivity such as that in iron pnicitides. 


PACS numbers: 74.25.Uv, 74.20.Rp, 74.70.Xa, 74.20.De 

I. INTRODUCTION 

ij/'J Vortex with In phase winding is a hallmark 
of macroscopic quantum state such as superfluidity and 
superconductivity^--. In superconductors, a vortex is accom¬ 
panied by a quantum of magnetic flux Oq = hcjle in a 
closed path with zero supercurrent. Since the quantization 
of magnetic flux is intimately related to the phase winding, 
superconductivity gap functions carrying intrinsic phase vari¬ 
ation induced by unconventional pairing symmetry should 
leave unique consequences on the response to external mag¬ 
netic held. Several interesting examples are available. A 
tricrystal ring of cuprate superconductor YBa2Cu307_(5 was 
observed to carry a half flux quanta Oo/2, which is the sig¬ 
nature for J-wave pairing symmetry^. In a ring-shaped setup 
composed of Nb and NdFeAsOo.88Fo.i2, flux jumps in odd- 
number multiple of half flux quanta were observed^-, which 
provide support to the S +_ pairing symmetry for iron-pnictide 
superconductors^—. Half-valued fluxoid jumps in magneti¬ 
zation curve of a thin annular coil composed of Sr2Ru04 were 
reported to be consistent with the / 7 -wave pairing symmetryii. 

The degree of freedom of relative phase difference in a two- 
component superconductor was first discussed by LeggetU^, 
and experimental observations on the collective Leggett mode 
were reported for MgB 24^*22. Phase solitons in one of the two 
superconducting gap functions associated with fractional flux¬ 
oid jumps have been investigated^"—. 

Another interesting possibility in multi-component super¬ 
conductors was raised some time ago 2 i whereby Josephson- 
like inter-component repulsions among three condensates 
can generate frustrations in phases of gap functions and in¬ 
duce a time-reversal-symmetry-broken(TRSB) superconduct¬ 
ing state even without external magnetic held, which corre¬ 
sponds to a high dimensional irreducible representation of 
point group similarly to the / 7 -wave superconductor—. Due 
to the discovery of iron-pnictide superconductors where sev¬ 
eral orbitals of Fe contribute to multi superconducting con¬ 
densates, this possibility becomes realistic and a considerable 



FIG. 1. Schematic setup with a loop of time-reversal-symmetry- 
broken (TRSB) superconductor, where the two halves are occupied 
by the two degenerate states with opposite chiralities. The three ar¬ 
rows denote the phases of order parameters and the yellow circle 
indicates the chirality referring to the mutual phases red ^ blue ^ 
purple among three condensates. Between the two halves of the loop 
there are two domain wall I and II accommodating inter-component 
phase kinks. 


amount of subsequent works have been devoted to discuss its 
thermodynamic stability and various novel properties'"^. 

Recently it has been revealed by the present authors 
that, in a Josephson junction between a conventional single¬ 
component superconductor and a multicomponent supercon¬ 
ductor in the TRSB state, the critical current should be asym¬ 
metric with respect to the current direction as the consequence 
of broken TRS^. As a matter of fact, unequal critical currents 
in opposite current directions were observed experimentally 
in a Josephson junction between Pbln and BaFei.8Coo.2As2^. 
Therefore, the TRSB state may have been realized already in 
iron-based superconductors, which is consistent with a micro¬ 
scopic analysis where band structures and strongly correlated 
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FIG. 2. Illustration of relation between the phase kinks at domain 
wall I and II and the flux trapped in the superconductor loop: (a) 
phase kinks between a same pair of components which can only 
trap integer multiples of flux quantum; (b) phase kinks between 
two different pairs of components which can trap a fractional flux 
(0 < 77 < 1). Colored curves at the domain walls refer to the phase 
variations in corresponding condensates, and D/y is for the gauge- 
invariant phase kink between condensate i and j. The dashed circle 
denotes the direction for counting phase winding in the loop. 
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FIG. 3. Schematic magnetization curve with fractional flux plateaus 
displayed together with gauge-invariant phase kinks denoted by 
[DikfDjk] in the order of domain walls I and II. 


effects are taken into account^. Cross checking this novel 
superconducting phenomenon becomes an important issue. 

In the present work, we address a new phase-sensitive prop¬ 
erty of the TRSB superconducting state. As schematically 
shown in Fig.[TJ we consider a loop of a multicomponent su¬ 
perconductor where the two halves are occupied by two TRSB 
states carrying opposite chiralities, accompanied by two do¬ 
main walls associated with inter-component phase kinks. We 
reveal explicitly that fractional flux plateaus appear in mag¬ 
netization curve associated with free-energy minima, where 
the domain walls accommodate phase kinks among different 
components leading to 271 phase winding along the loop only 
in one or two of the three condensates. While the heights 
of fractional flux plateaus depend on material parameters and 
temperature, they form pairs with heights related by the flux 
quantum Oq, which is a unique signature of the TRSB super¬ 
conducting state and can be used to confirm the state itself. 
In a more general point of view this provides a novel chance 
to explore relative phase difference, phase kink and soliton in 
ubiquitous multicomponent superconductivity. 

The remaining part of this paper is organized as follows. 
The mechanism for fractional flux plateaus in presence of do¬ 
main walls is discussed in Sec. II. We then simulate the mag¬ 
netizing process with the time-dependent Ginzburg-Landau 
(TDGL) approach and confirm fractional flux plateaus associ¬ 
ated with free-energy minima in Sec. III. In Sec. IV, we study 
an asymmetric loop where the widths of two domain walls 
differ. In Sec. V we discuss the stability of the domain-wall 
structure and temperature dependence of the fractional flux 
plateaus. Finally we give a summary in Sec. VI. 


II. FRACTIONAL FLUX IN SUPERCONDUCTING LOOP 

In order to reveal the essence of physics we first con¬ 
sider an isotropic TRSB state which is generated by three 


equivalent condensates with equal mutual repulsion. For 
simplicity all order parameters are taken as v-wave from 
now on. The two states in the loop are given by 'F = 
1/73} = ^«4;r/3| opposite chirali- 

ties (see Fig.[T]). Across each of the two domain walls between 
the left and right halves of the superconductor loop, there is 
a phase kink where the intercomponent phase difference be¬ 
tween two of the three order parameters shrinks to zero and 
reopens in the opposite way continuously, resulting in a sign 
reversal in the phase difference at the two sides of domain 
(see Fig. 2 ). We notice that phase kinks are gauge-invariant 
objects, which inevitably appear at the interface between two 
bulks of TRSB states with opposite chiralities. 

When the two domain walls accommodate the phase kink 
between the same two condensates, such as that between con¬ 
densate 1 and 2 defined as D12 in Fig.^a), the phase rotation 
integrated in a counterclockwise manner (indicated by L in 
Fig.O over the two domain walls cancel each other, resulting 
in the same phase winding in all the three condensates. In this 
case, the flux trapped in the loop is an integer multiple of flux 
quantum Oq when the loop is thick enough to fully screen the 
magnetic field. 

The situation differs when the two domain walls accommo¬ 
date different phase kinks, such as D12 and D23 in domain wall 
I and II respectively shown in Fig. Ob). By inspection one 
sees that if/2 rotates Ttt/S anticlockwise over the two domain 
walls, while and 1I/3 rotate - 27 r/ 3 . When the external mag¬ 
netic field provides the additional phase rotation of It: 12 in all 
condensates, a state with It: phase winding in if/2 and 0 in both 
ij/i and j/^3 is stabilized. This yields a fractional flux quanta 
O0/3 in the loop. The state with a fractional flux trapped in 
this loop is expected to be stable in a certain range of exter¬ 
nal magnetic field, which leads to a fractional flux plateau in 
magnetization curve shown schematically in Fig.O 

The above discussion can be elucidated by the integration 
of magnetic flux over the superconductor loop using the GL 
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formalism where the supercurrent is given 



( 1 ) 


with ruj and (fj the effective mass and phase of component- 
j. For a thick loop, the supercurrent is zero deep inside the 
superconductor. In this case the magnetic flux trapped in the 
loop is given by the line integration of phase differences as 
can be seen from Eq. ([T]) 


O = 

271 


Iti 


i 


PlV(pi + + / 73 V ^3 


dl 


(f)V(pidl+ j 
Jc Jdw 


/7l + /72 + P3 

PlS<P\l + P3^<P\3 
/7l + P2 + P3 


dl 


( 2 ) 


with pj = \il/j\^lmj and (fij = (fj - (pi for /, 7 = 1,2, 3, where 
is a closed path along the loop with zero supercurrent 
everywhere (see Fig. 1), and the ”DW” denotes domain-wall 
regimes (gray parts in Fig. 2) with phase kinks. In the second 
line, we divide the integrand into two terms, indicating two 
contributions to the total magnetic flux. The first contribution 
should be an integer multiple of It: due to the single-valued 
wave function in the loop. The integrand in the second contri¬ 
bution is nonzero only on domain walls. This contribution is 
nonzero when two different phase kinks are realized at domain 
wall I and II, with the value depending also on the quantities 
Pj’ 

Presuming the same length of domain walls I and II, the 
two configurations [Dtk/Djk] and [Djk/Dik] at the domain walls 
[I/II] take the same free energy. However, integrating phase 
differences for these two configurations along the closed path 
in the counterclockwise manner (see Fig. O results in op¬ 
posite fractional values of 2;r in the second term in Eq. 
Therefore, these two configurations give two fractional fluxes 
Oi and O 2 related by the flux quantum Oq. Fractional flux 
plateaus with corresponding configurations of phase kinks are 
schematically shown in Fig.O 


III. TDGL APPROACH 


Here we adopt the GL formalism to check the thermody¬ 
namic stability of states carrying fractional fluxes. The GL 
free-energy functional of a three-band superconductor with 
Josephson-like inter-component couplings is given by^^*^ 


F= Z 

7 - 1 , 2,3 




Z yjk + C.C.\ + i (V X A)^ , (3) 


where aj is a temperature-dependent coefficient which is neg¬ 
ative when T < Tcj and positive when T > Tcj, with Tcj 
the critical point of the superconducting component- j before 
considering intercomponent couplings, and jjk is an inter¬ 
component coupling taken as constant for simplicity. For 
T 12 T 13 T 23 < 0, a TRSB superconducting state appears when 




FIG. 4. Amplitudes and phases of order parameters in two TRSB 
states with parameters a\ - 0 . 012 , 0^2 = 0.013, = 0.011, 7 J 2 = 

723 “ -0.24, 7'3 = -0.25, m\ - - I, - 1.1, - 

P' 2 , - ^ and K\ - 1.5. See text for definitions of the dimensionless GL 
parameters. 


the coefficients in Eq. ([3]) satisfy conditions revealed in a pre¬ 
vious work^. To be specific, we put 712 , 713,723 < 0, namely 
all repulsive Josephson-like couplings, since it is easy to see 
that other TRSB states can be generated from this case by a 
simple gauge transformation. 

We adopt dimensionless quantities given by^ 

X = ^i(0)x', A = ii(0)H,,i(0) V2A', J = 

osj = 7jk = -O'i(0)r}i, 

nij = mim/, ki = /li(0)/^i(0), F = GqF' 


with t^io(O) = -ai(0)lj3u ^i(O) = ^niic^/ [47r^^Q(0)(2e)2], 

^i(O) = ^J-n^/[2mia\(0)], Htc\ = and 

Go = (0)/47r. In the dimensionless units, the GL free 

energy is rewritten as 


F'= Y. 

7 - 1 , 2,3 



V 

<1 

1 

> 

|2l 


\\lKl 7 ' 

^1 


J,k= 


,2,3\j<k 


( 4 ) 


The system can be described by the following TDGL equa¬ 
tions in the zero-electric potential gauged 


dll/'. 1/1 

^ = -v-A'U;+ ^ 

’ k=\, 2 Xk*i 

( 5 ) 

with j = 1,2,3 and 

Z Vv);.-A')-VxVxA' (6) 


IKi 


with cr the coefficient of normal conductivity. At equilibrium 
the left-hand sides of Eqs. ([5]) and ^ are zero, which gives 
the GL equations. By solving three GL equations in Eq. ^ 
with A' = 0, we obtain the amplitudes and phases of the con¬ 
densates at zero magnetic field^. 

In the present multi-component superconducting system, 
the critical point Tc is higher than Tcj for j = 1,2,3^. In 
order to make sure that the GL formalism is valid for investi¬ 
gating the thermodynamics properties of the coupled system, 
we choose temperature satisfying Tcj < T < T^ where a'j 
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H[Htci(0)] 


FIG. 5. (a) Magnetization curve of thermodynamically stable state 
and (b) Gibbs free energies of several competing states in the su¬ 
perconductor loop shown in Fig. 1 upon sweeping external magnetic 
field. Parameters are the same as Fig. |4| The superconductor loop 
is of square shape with outer frame of 24/Ii(0) and inner frame of 
8Ti(0). 


are positive and small. In Fig. |4] we display the two TRSB 
states with opposite chiralities, which are used for the follow¬ 
ing study on the magnetization process. 

We take a square loop with external dimension 24/li(0) x 
24/li(0) and the wall thickness 8/li(0) to investigate the mag¬ 
netic response. From Eq. ([T]) we obtain the the penetration 
length in the dimensionless form. 

For the states given in Fig.|4l A = 1.22/li(0) is much smaller 
than the thickness of loop wall. Therefore, one can take the 
closed path ”C” as the middle line of the superconductor loop 
where supercurrent is negligibly small (see the Appendix). At 
the edge of superconductor, we presume that no supercurrent 
flows out of the superconductor, and the B fleld at the external 
edge of superconductor loop is flxed to the value of applied 
magnetic fleld. 

The magnetization curve derived from TDGL equations ([5]) 
and (O is shown in Fig.[5l where fractional flux plateaus cor¬ 
responding to states with free-energy minima are obtained. 
When the magnetic fleld is small, there is no flux pene¬ 
trating into the loop. As the magnetic fleld increases to 
HIHtcxijA) = 0.009, the stable state takes di fractional flux 
Oi = 0.260o. This state remains stable until the magnetic 
fleld = 0.019, yielding a plateau in the magneti¬ 

zation curve. As the magnetic fleld increases further, another 
state appears with fractional flux O 2 = 0.740o in the regime 
0.019 < HIHtciif)) < 0.029. We can see that both Oi and O 2 



FIG. 6. (color online) Magnetization curve with two fractional flux 
plateaus for an asymmetric superconductor loop, where the width of 
side including domain II is enlarged to l2Ai(0) from that given in 
Fig. [ 5 ] The parameters are the same as Fig.|4]except for a[ = 0.025, 
= 0.028 and a'^ = 0.022. 


deviate from Oo/3 and 20o/3, because of the three inequiv¬ 
alent condensates in the system. Nevertheless, it is clear that 
the relation Oi -r O 2 = Oq is satisfied as revealed in Sec. II. For 
even larger magnetic fields, the stable state permits one flux 
quantum Oq inside the loop. Note that the magnetic fields for 
the fractional flux plateaus are very small as compared with 
the typical field Htd (0) and thus there is no vortex inside the 
body of the superconductor. 

We check the phase kinks at the two domain walls I and II 
and the phase windings along the loop in the three conden¬ 
sates at the fractional flux plateaus. At Oi = 0.260o, phase 
kinks Di 2 and D 23 are realized at regions I and II respectively, 
and 0^2 rotates In along the loop leaving 0i and 03 unwinding. 
At Oi = O. 74 O 0 , the phase kinks D 12 and 7)23 are at domain 
walls II and I, in contrast to the case of Oi = 0.260o, and 0i 
and 03 rotate In with 02 unwinding. At integer flux quanta 
0 = 0 and O = Oq, the phase kink is D \2 at both domain 
walls I and II. All these are in accordance with the discus¬ 
sion in Sec. II. In general, there are at most six fractional flux 
plateaus between the integer flux quanta 0 and Oq. For the 
GF parameters given in Fig.|4]we can only see two plateaus 
in Fig. Oa) because the free energies of domain walls satisfy 
F{JA\ 2 ) ^ F{JA 2 f) < F{D\ 2 ,) such that only the phase-kink pair 
Di 2 and D 23 is stabilized. 


IV. ASYMMETRIC SUPERCONDUCTOR LOOP 

Up to this point, the superconductor loop is presumed to 
have the same width at domain wall I and II, which guaran¬ 
tees the degeneracy of domain-wall energy between phase- 
kink pairs [Dik/Djk] and [Djk/JAik ]. In this case, it is easy to 
see that the magnetization curve in Fig. [3] should be symmet¬ 
ric with respect to the direction of magnetic field. In general, 
the two widths can be different. In the latter case [DiflDjk] 
may be unstable even though [Djk/Dik] is stable and associ¬ 
ated with the free-energy minimum. As shown in Fig.[6]for an 
asymmetric superconductor loop, the fractional flux plateau at 
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FIG. 7. (color online) Gibbs free energy for the system in Fig. 1 with 
two halves occupied by the two TRSB states with opposite chiralities 
as a function of the location of domain walls x defined in the inset. 
Parameters and sample size are the same as Fig. 5. 


O = 0.740o remains stable for 0.017 < HjHtcxi^) < 0.023, 
while that at O = 0.260o disappears in contrast to Fig.[5l since 
they are associated with different phase-kink configurations at 
domain wall I and II. It is worth noticing that even in this 
asymmetric loop the plateau at O = -0.260o is still stable for 
-0.020 < < -0.014, since it is associated with the 

same phase-kink configuration with that at O = 0.740o and 
a difference of flux quantum Oq coming from the first term in 
Eq. (|2l). In this asymmetric loop, the magnetization curve is 
asymmetric with respect to the direction of the magnetic field 
as in Fig.[6l The property that fractional flux plateaus in posi¬ 
tive and negative magnetic fields are paired with the difference 
of flux quantum Oq is robust, and can be taken as a crosscheck 
for fractional flux plateaus originated from the TRSB state. 


V. DISCUSSIONS 

In the present work we study the case that the left and right 
halves of the superconductor loop take the two TRSB states 
with opposite chiralities. This situation can be realized in 
experiments by cooling the whole system from temperature 
above with laser heat pulse irradiated on regions I and II^. 
The two halves condensate independently and by chance ar¬ 
rive at the different TRSB superconducting states, leading to 
the two domain walls at region I and II after releasing the ir¬ 
radiation. In order to check the stability of this configura¬ 
tion, we estimate the free energy of the whole system in terms 
of the TDGL approach. As shown in Fig. [71 the state with 
two domain walls located at the middle of the top and bottom 
sides of the loop corresponds to a free-energy minimum. The 
domain walls once generated should be stable since moving 
them outside the loop is prohibited by a large free-energy bar¬ 
rier which is produced by an elongated, single domain wall 
during the process of domain-wall relocation (see that at X 2 in 
the inset of Fig. [7]). The stability of the present setup against 
relocating one of the two domain walls along the loop can be 


provided by widening the left and right arms of the loop. The 
increase in free energy in states with fractional fluxes and inte¬ 
ger flux quanta upon application of external magnetic field is 
smaller by one order of magnitude than the free-energy barrier 
as seen in Fig. [71 which justifies the discussion on fractional 
flux plateaus in the present work. 

At this point we notice that the free-energy barrier in Fig. 7 
generated by the two TRSB states at the two halves of su¬ 
perconductor loop is crucially important for the thermody¬ 
namic stability of fractional flux plateaus. From Eq. (O one 
might think that a two-component superconductor or a three- 
component one with preserved TRS can also accommodate 
fractional fluxes. However, in these cases there is no free- 
energy barrier like that in Fig. 7, and states with fractional 
fluxes are unstable. 

The amplitudes and inter-component phase differences of 
order parameters change with temperature, leading to varia¬ 
tion in stable domain-wall structures. As a result, both height 
and width of a fractional flux plateau should change as tem¬ 
perature is swept. During field-cooling or field-heating pro¬ 
cesses, jumps between states with fractional fluxes and/or in¬ 
teger multiples of flux quantum can take place. 

Half-valued fluxoid jumps in magnetization curve were re¬ 
ported in a small multicomponent superconducting sample 
comparable with penetration lengtbi^. This results in incom¬ 
plete screening of magnetic field, and thus the magnetization 
curve exhibits a finite slope for any external magnetic field. 
Namely, there is no fractional flux plateau in their system. In 
previous studies on vortex states of TRSB superconductor it 
was discussed that vortex cores of different condensates can 
deviate from each other in spac o^^i^^" — . However, without 
a closed path along which supercurrent is zero everywhere, 
there is no fractional flux plateau in the magnetization curve. 


VI. CONCLUSION 

To summarize, we have studied the magnetic response of a 
loop of three-component superconductor with two degenerate 
time-reversal symmetry broken states at two halves. When 
the two domain walls between the two halves accommodate 
different intercomponent phase kinks, fractional flux plateaus 
appear in the magnetization curve which form pairs related to 
each other by the flux quantum. These properties are expected 
to be helpful for detecting experimentally the time-reversal 
symmetry broken superconducting state which can be realized 
in iron-pnictide superconductors. In general, this endeavour 
provides a novel chance to explore relative phase difference, 
phase kink and soliton in ubiquitous multi-component super¬ 
conductivity. 
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FIG. Al. (color online) Distribution of supercurrent density J' in the 
superconductor loop for the fractional flux plateau of O = 0.260o in 
Fig. 5, and that of order parameters on domain wall II. The directions 
and lengths of arrows refer to phases and amplitudes of order param¬ 
eters, where |j^ 2 l and are suppressed to 0.63 and 0.77 of the bulk 
values while |j/^i | remains almost unchanged. 
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APPENDIX: DISTRIBUTION OF SUPERCURRENT IN THE 
SUPERCONDUCTOR LOOP 

An example of the distribution of supercurrent density in 
the superconductor loop and order parameters on domain 
walls is shown in Fig. lAIi where the superconductivity sur¬ 
vives in all components and the Meissner effect screens the 
magnetic field and thus the total supercurrent to zero. The 
important feature here is that, although the supercurrents in 
individual components are not zero due to the phase shifts, 
they flow in opposite directions and cancel each other, lead¬ 
ing to zero total supercurrent. The situation of a domain wall 
as discussed in the present work differs considerably from a 
vortex, where the phase winding induces a divergent kinetic 
energy, which necessitates total suppression of the amplitude 
of the superconducting order parameter at the vortex core. 
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